Numerical study of persistence in models with absorbing states 



O 

o 
o 

(N 
•*-> 
O 

O 



(N 
O 
O 

o 
o 

-I— > 



c 

o 
o 



X 



Ezequiel V. Albano 1 and Miguel A. Munoz 2 
Instituto de Investigaciones Fisicoquimicas Tedricas y Aplicadas, (INIFTA), CONICET , UNLP, GIG (Bs.As.). Sucursal 4, 

Casilla de Correo 16, (1900) La Plata, ARGENTINA. 
2 Institute Carlos I for Theoretical and Computational Physics 
and Departamento de E. y Fisica de la Materia, 
Universtdad de Granada 18071 Granada, Spain. 
(February 1, 2008) 

Extensive Monte Carlo simulations are performed in order to evaluate both the local (#;) and 
global (9 g ) persistence exponents in the Ziff-Gulari-Barshad (ZGB) ( Phy. Rev. Lett. 56, 2553, 
(1986) ) irreversible reaction model. In two dimensions and for the second-order irreversible phase 
transition (IPT) we find, that both the local and the global persistence exhibit power-law behavior 
with a crossover between two different time regimes. On the other hand, at the first-order IPT, 
characteristic of the ZGB, active sites are short lived and the persistence decays more abruptly, not 
being clear whether it shows power law behavior or not. In order to analyze universality issues, we 
have also studied another model with absorbing states, the contact process, and evaluated the local 
persistence exponent in dimensions from 1 to 4. A striking apparent super-universality is reported: 
the local persistence exponent seems to coincide in both one and two dimensional systems. 



I. INTRODUCTION 

An avalanche of articles devoted to the study of persis- 
tence in different physical contexts has appeared recently 
in the literature. From a theoretical point of view, the 
appearance of new nontrivial critical exponents as: i) the 
"short-time" exponent, A, needed to characterize the two- 
time correlation functions in systems relaxing in the pro- 
cess of quenching from infinitely high temperatures to the 
critical temperature T c [jjjpl , and ii) the global persistence 
exponent 8 g , related to the probability P g (t) oc t~ 0g , that 
the global order parameter, (e.g. the magnetization in a 
ferromagnet), has not changed sign up to a certain time 
t after a quench to T c S, constituted an exciting sur- 
prise. A number of theoretical, computational, and even 
experimental analysis have been performed afterwards. 
From the theoretical viewpoint it has been shown that 
for processes where the global order parameter dynamics 
is Markovian, that the global persistence exponent can 
be related to other "traditional" critical exponents: in 
particular, the following scaling law holds 
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non-Markovian processes (perturbing around a Gaussian 
and Markovian process |4||§]). 

Recently, other similar non-trivial exponents have been 
undercovered and measured in different systems. Some 
of them are: 

• (i) The local persistence probability Pi oc t~° l de- 
fined as the probability that the local order param- 
eter at a given point x has never changed sign after 
the initial time (e.g. the probability than one spin 
has never been flipped in a ferromagnet) H^). Lo- 
cal persistence exponents have been measured in 
real systems as liquid crystals or soap bubbles (see 
jl| and references therein). 

• (ii) The block persistence probability Pl oc t~ 0L 
defined as the probability that the order param- 
eter integrated over a block of linear size L has 
never changed sign since the initial time in a phase- 
ordering process at finite temperature (l(],[ll ] ; 

• (iii) The generalized persistence probabilities and 
exponents introduced by Dornic and Godreche 
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where d is the dimensionality and r\ is the static criti- 
cal exponent of the order parameter correlation function. 
Nevertheless, the required hypothesis is not typically sat- 
isfied for most of the usually studied statistical models; 
i.e. the dynamics of the global order parameter can be 
argued to be gencrically non Markovian. In these more 
generic cases Eq. (1) does not hold and g is a new inde- 
pendent non-trivial exponent j3j. In general, persistence 
exponents depend on the system evolution as a whole 
and, therefore, analytical estimations of them are scarce 
and difficult [gj-Q- O ne important theoretical contribu- 
tion is by Majumdar and Sire, who proposed a method 
to calculate autocorrelation functions perturbatively for 



• (iv) The persistence exponents for domains. Found 
first by Krapivsky and Ben-Nairn in Ising system 
and sub seq uently in diffusion by Majumdar and 
Cornell 13fl. These can also be generalized to other 
different types of pattern: each pattern decays with 
a different exponent, a single spin or the global spin 
being just two specific patterns. 

• (v) The persistence with partial survival first stud- 
ied by Majumdar and Bray and subsequently by 
other authors jlj]. Recently it has generated a lot 
of interests in the context of inelastic collapse in 
granular materials. 
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For a review of persistence studies in reversible systems 
see e.g. OgQ. 

Even though considerable effort has been devoted to 
the study of persistence in various models exhibiting re- 
versible phase transitions (see Jl|-||||[)] [ll 16 2(]] and ref- 
erences therein) as well as in diffusion-reaction systems 
(see 21 and references therein), far less numerous 
are the studies of persistence in systems exhibiting ir- 
reversible phase transitions (IPT's) and, in particular, 
in systems presenting a critical point separating an ac- 
tive from an absorbing phase (some exceptions can be 
found in j23-E9|}. In this context, recently, Hinrichsen 
and Kodusely |23|] have performed a numerical study of 
the local and global persistence in (1 + 1)— dimensional 
directed percolation (DP) |^6| (see also ^7j). They found 
that the local persistence probability at the critical point 
(separating the absorbing from the active phase) decays 
as a power law with an exponent 9i DP ~ 1.5. Also, global 
persistence measurements seemed to be consistent with 
9 g DP ~ 1.5, a result that has to be taken with some 
caution since, on the one hand, the evaluation of 9 g DP 
requires extensive simulations (see the criticisms to this 
result in [p4| ) and on the other hand, in all the previously 
known models (mainly for equilibrium critical systems) 
g < 6i j3). Furthermore, it was conjectured that both 
9 g and 61 are indeed universal exponents characterizing 
the DP universality class. Another interesting point, is 
that the persistence exponents in this case seem to be 
independent on the initial condition at odds with what 
happens in other well known cases as, for example, the 
two-dimensional Ising model with Glauber dynamics Q . 
Finally, Hinrichsen and Koduvely succeeded in relating 
persistence exponents in directed percolation to certain 
"return" probabilities with an absorbing boundary or an 
active source p0[ . Another important contribution to 
the study of this IPT was presented in the recent pa- 
per by Oerding and van Wijland, in which the global 
persistence exponent has been calculated analytically by 
combining the perturbative method developed by Ma- 
jumdar and Sire ^] with standard renormalization group 
techniques p4fl . These same authors have thrown some 
doubts on the validity of the numerical estimations ob- 
tained by Hinrichsen and Koduvely (as will be discussed 
in section III). 

In order to contribute further to the understanding 
of persistence in systems with irreversible phase transi- 
tions and try to shed some light on the aforementioned 
issues, the present manuscript is devoted to the numerical 
study and evaluation of persistence exponents of the Ziff- 
Gulari-Barshad (ZGB) model for the catalytic oxidation 
of CO 1 28 1 , as well as in the contact process 29 , 2(| in var- 
ious dimensions. The ZGB model in d = 2— dimensions 
has a twofold advantage: (i) it exhibits a second-order 
IPT unambiguously placed at the DP universality class 
p0| ], and (ii) it also exhibits a first-order IPT |^8| where, 
as expected, the system does not show scaling neither 
universal behavior M. Therefore, the ZGB model pro- 



vides a suitable framework for the study of the persis- 
tence in both, first and second order irreversible critical 
points. Furthermore, the present study extends the in- 
vestigation of persistence exponents to higher dimensions 
in IPT's, by performing numerical studies of the contact 
process. 

The manuscript is organized as follows: In section II 
the ZGB model and the simulation method are described 
in detail. In section III we define the local and global 
persistence, while section IV is devoted to the presenta- 
tion and discussion of the main results for both the ZGB 
model and the contact process. Finally in section V the 
main conclusions are presented. 



II. THE ZGB MODEL AND THE SIMULATION 
METHOD 

The ZGB model mimics the catalytic oxidation of car- 
bon monoxide on a transition metal surface ]28|] , namely 
2CO + 2 -»■ 2C0 2 . The lattice-gas version of the ZGB 
model is also known as the monomer-dimer model (MD). 
where A = CO is the monomer because it needs a sin- 
gle adsorption site on the surface, and B 2 = 2 is the 
dimer which adsorbs dissociatively, and consequently re- 
quires two neighboring sites on the catalyst surface to 
be adsorbed. It is assumed that the monomer-dimer re- 
action proceeds according to the Langmuir-Hinshelwood 
mechanism: 



A(g) + S^A(a), 



B 2 {g) + 2S^ 2B(a), 



A(a) + B(a) -> AB{g) + 25, 



(2) 



(3) 



(4) 



where 5 is an empty site on the surface, while (a) and 
(g) refer to the adsorbed and gas phase, respectively. 

The MD model uses a square lattice to represent the 
catalytic surface. The Monte Carlo algorithm for its sim- 
ulation is as follows : i) an A or B 2 molecule is selected 
randomly with relative probabilities Ya and Yb, respec- 
tively. These probabilities are the relative impingement 
rates of both species, which are proportional to their par- 
tial pressures. Due to the normalization, Ya + Yb = E 
the model is characterized by a single parameter, say Ya- 
If the selected species is A, a surface site is selected at 
random, and if that site is vacant, A is adsorbed on it 
[Eq.(2)]. Otherwise, if that site is occupied, the trial ends 
and a new molecule is selected. If the selected species is 
B 2 , a pair of nearest-neighbor sites is randomly chosen 
and the molecule is adsorbed on them only if they are 
both vacant [Eq.(3)]. ii) After each adsorption event, 
the nearest neighbors of the added molecule are exam- 
ined in order to account for the reaction given by Eq. (4) . 
If more than one pair [B(a), A(a)] are identified, one of 
them is randomly selected and removed from the surface 
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(for more details on the MD and the simulation technique 
see, e.g. HH). 

Interest in the MD model arises due to its rich and 
complex phenomenology. In fact, in two dimensions and 
for the asymptotic regime (t — * oo), the system reaches a 
stationary state whose nature depends solely on the pa- 
rameter Y A . In fact, for Y A < Y 1A = 0.387368 (Y A > Y 2A 
= 0.52554) the surface becomes irreversibly poisoned by 
B (A) particles. On the other hand, for Y\ A < Y A < Y%a 
a steady state with sustained production of AB is gener- 
ated. Figure 1 shows plots of the rate of AB production 
(R A b) and the surface coverages of A (0 A ) and B (9b) 
versus Y A , respectively. Just at Y 1A and Y 2A the MD 
model exhibits irreversible phase transitions (IPT's) be- 
tween the reactive regime and poisoned states, which are 
of second and first order respectively. The second-order 
IPT belongs to the universality class of directed percola- 
tion p6| and is rather well understood ]30[ ] . Furthermore, 
as it is shown in Figure 1, when Y A increases towards Y 2A 
the catalytic activity increases, but when Y 2A is reached 
large A clusters suddenly emerge and cover the whole 
lattice. The transition occurs abruptly, with discontinu- 
ities of the coverages and activity, unveiling its first-order 
nature. 



III. LOCAL AND GLOBAL PERSISTENCE 

In a system, such as directed percolation, possessing 
absorbing states, the local persistence probability can 
be measured by starting the system with a homogeneous 
random initial conditions, and evaluating the probabil- 
ity that a given, originally inactive, site has not become 
active up to time t pOj ]. It should be noticed that this 
quantity depends upon the whole system evolution up to 
time t, and therefore it can be enviewed as an infinite- 
point correlation function in the time direction. Conse- 
quently Pi (t) is a rather non-trivial quantity, difficult to 
study analytically as discussed before. Numerical sim- 
ulations in a (1 + l)-dimensional DP model shows that 
Pi (t) decays algebraically at criticality: 



Pi(t) oc r 



(5) 



with B{ w 1.5 p6| . In the case of the ZGB model, and for 
the continuous phase transition, the absorbing state cor- 
responds to the surface of the catalyst fully covered by 
B-species. Such species plays the role of inactive sites 
(for the other transition, i.e. the first-order one, the 
roles of the A and B particles are exchanged). We de- 
fine P^ GB (t) as the probability that a given inactive site 
does not become active up to time t. Consequently, for 
the second-order (first order) transition simulations are 
started with lattices partially covered at random by A- 
species (-B-species) (9 A (B)(t = 0)). The only mechanism 
capable to activating a site is the reaction A + B — > AB 
as specified by equation 4. 



We have numerically evaluated the persistence distri- 
bution function Df GB (t), that is, the probability for a 
given site to become active in the time interval between 
t and t + dt. Accordingly, one has 
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ZGB 



(t) (X t 
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J + l) 



(6) 



where 9f GB is the persistence exponent (Di ZGB (t) is 
simply the time derivative of P ZGB ). 

In systems exhibiting reversible phase transitions, the 
global persistence is usually defined as the probability 
that the global order parameter (e.g. the total magne- 
tization in the Ising model) does not change sign up to 
time t. In contrast, for systems with absorbing states, 
such as DP, the global order parameter given by the 
density of active sites p A s(t), is a strictly positive quan- 
tity and therefore the standard definition is not applica- 
ble. Hinrichsen and Koduvely |20| have proposed that, 
instead, one can evaluate the probability that the de- 
viation of the order parameter from its mean value, 
&pAs(t) = PAs(t)~ < PAs{t) >, does not change sign up 
to time t. In the present study, and close to the second- 
order IPT, we have considered 



PAs(t) 



1 



(7) 



where 9y, 9 A and 9b are the density of empty and occu- 
pied by A and B— species sites, respectively. However, in 
finite systems as directed percolation the deviation prob- 
ability depends on the sign of Ap A s(t). This asymmetry 
should vanish in the thermodynamic limit where Ap A $(t) 
becomes a Gaussian process. In finite systems, however, 
the renormalized variance of the fluctuations is no longer 
constant but increases with the value of Ap A s(t) causing 
the effective time scales for positive and negative fluctua- 
tions to be different ||||. Therefore, P ( g ~ ] (t) {Pg +) (*)) 
is defined as the probability of Ap A s(t) to remain nega- 
tive (positive) from an initial time U n up to time t. We 
have also evaluated the distribution function Dg GB (t), 
i.e., the probability that Ap A s(t) changes sign in the 
time interval between t and t + dt. In order to evaluate 
Ap A s(t) we have first to calculate a "calibration curve" 
given by p A s oc t~ 9 where 9 ~ 0.4505 is the correspond- 
ing DP exponent in (2+1)— dimensions [ p3p^ ], describing 
the time decay of a homogeneous initial condition at crit- 
icality: every-time p A s(t) intersects this curve, Ap A s(t) 
changes sign. 

A second important issue related to the global persis- 
tence exponent, worth to discuss before proceeding, is the 
time regime in which it should be measured. Following 
pij , the global exponent is well defined in the regime in 
which the dependence on the initial state has been erased 
and, on the other hand, there is an upper bound to the 
validity of measurements, given by the upper cutoff in- 
duced by the finite system size. Putting together these 
two constraints |24f| 
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(for general systems with absorbing states, one has to 
substitute 8b by the initial density of active sites). 

Using the DP values for the exponents, this implies 
18300 > t > 1 for system size L = 256, and 65000 > 
t 1 for L = 512. Therefore, as a safe, conservative 
limit, we start measuring at tj n = 100. 

It is worth stressing that in order to obtain reliable 
data for persistence exponents one has to perform exten- 
sive simulations. The evaluation of the local persistence 
requires long ranges (up to 10 4 Monte Carlo steps), and 
reliable statistic is obtained averaging over 10 4 (2.5 xlO 3 ) 
different initial configurations for L = 128 (L = 256). 
On the other hand, each single point for the global per- 
sistence distribution requires a whole run; therefore good 
statistics is much tougher to obtain that for local persis- 
tence. We averaged over 10 6 and 2.10x 5 different runs 
for L = 256 and L = 512, respectively. 



IV. RESULTS AND DISCUSSION. 

In this section we report our main numerical findings. 
The value of the second order critical point of the ZGB 
model reported previously, Ya = 0.387368 |50|, corre- 
sponds to the infinite size limit; while for the system 
sizes we will consider, namely 128 * 128 and 256 * 256, 
finite size corrections shift the critical value to slightly 
larger values: in particular, all the simulations at the 
critical point reported in what follows are performed at 
Y A = 0.3907 f(|. 

Figure 2 shows a log-log plot of Di(t) versus t at 
the second-order IPT. Two regimes can clearly we ob- 
served: i) The short time one, for t < 200, where Di(i) 
exhibits a power law behavior Di(t) oc t~( 9l+1 \ with 
0i ~ 1.00 ± 0.05, and (ii) the asymptotic regime, for 
t > 200 , with the exponent 6f GB ~ 1.50 ± 0.10. Since 
the phase transition belongs in the DP universality class 
p0| , it is quite surprising that the obtained exponent, 
within error bars, is almost the same than the one re- 
ported by Hinrichsen et al. pG], for a DP process in 



(1 + 1)— dimensions, namely 9f 



1.50 ± 0.03. This 



finding suggests that for DP, the difference between the 
persistence exponents in (1 + 1)— and (2 + 1)— dimensions 
is very small. In addition, we have also verified the lack of 
dependence of these exponents upon modifications of the 
initial condition, as shown in figure 3. In fact, runs per- 
formed using 0.05 < Osit — 0) < 0.20 show clearly that 
neither the short-time regime nor the persistent regime 
depend on the initial density of inactive sites. This find- 
ing is in accordance with Hinrichsen et al. observation 

In order to gain some insight on the origin of the 
crossover observed in figures 2 and 3, we have evalu- 
ated the time dependence of the fraction of persistent 
sites with NN = 0, 1, 2, 3 and 4 empty nearest-neighbor 
respectively, as a function of time, as it is shown in fig- 
ure 4. It is clear that different NN values vanish at 



distinct characteristic times, e.g. t ~ 10, (NN = 0); 
t ~ 15, (NN = 1); t ~ 30, (NN = 2); t ~ 120, (NN = 3); 
and t ~ 300, (NN — 4), respectively. Therefore compar- 
ing the results of figures 2-3 and 4 we conclude that the 
crossover in figures 2-3 takes place when almost all persis- 
tent sites are blocked by inactive (absorbing) ones. The 
physical interpretation of the observed crossover is there- 
fore the following: There is an initial regime in which the 
persistence of a given site is dominated by its open neigh- 
boring sites: particles can land directly in the neighbor- 
ing sites causing desorption, and consequently the death 
of persistent sites. However, as the number of such open 
sites is reduced on average to a value close to zero, there 
is a crossover to the true asymptotic regime in which one 
has to wait for local rearrangements in order to reach 
persistent sites. 

In order to further explore the dependence of the per- 
sistence exponents on the dimensionality we have also 
performed extensive simulations f34jj , of another system 
with absorbing states belonging to the DP universality 
class, namely the contact process (CP) [^6 29 . In this 
model, proceeding analogously, we observe no crossover 
analogous to the one described for the ZGB, asymptotic 
regimes can be obtained with less computational effort, 
and the exponent values arc far more accurate. We simu- 
late the CP in dimensions from 1 to 4. The largest system 
sizes considered are 10 4 in d = 1, 256*256 in d — 2, 50 3 in 
d = 3, and 32 4 in d = 4. Our results are 9 t ~ 1.50 ± 0.01 
for (1 + 1)-, and 6> ; ~ 1.50 + 0.01 for (2 + 1)- dimensions 
Q (see Fig. 5 ). In fact, the two curves in Figure 5 are 
strikingly parallel. This fact confirms our rather strange 
finding, that the local persistence exponent is (almost?) 
the same in one and two-dimensional systems with ab- 
sorbing states. It could be the case, that there is some 
small difference between the exponent values in d = 1 
and d = 2, but within our numerical accuracy they are 
absolutely indistinguishable. The apparent coincidence 
between the results in one and two dimensions is quite 
intriguing, and so far, we do not have any satisfactory 
explanation for this surprising result. 

For the sake of completeness, and in order to test 
whether this apparent "super-universality" extends also 
to higher dimensions, we have also measured the local 
persistence exponents in d — 3 and d = 4 as said before. 
The results for 6» z are 1.33 + 0.03 in d = 3, and 1.15 + 0.05 
in d = 4, indicate a monotonous diminution of the local 
persistence as a function of dimensionality Q . 

It is well known that for equilibrium models the persis- 
tence exhibits generically a power-law behavior below the 
critical temperature |ll| , |l2| . In order to establish analo- 
gies with irreversible systems having absorbing states, we 
have also measured the local persistence away from crit- 
icality, i.e. within the reactive regime of the ZGB model, 
as it is shown in figure 6. For Ya = 0.40 (i.e. very close 
to the critical threshold) the power-law characteristic of 
the short-time regime is still observed, but the persis- 
tent regime shows a clear departure from the behavior 
observed at criticality. Going deeper into the reactive 
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regime, e.g. for Ya — 0.45, the plot shows a marked cur- 
vature and the existence of a power-law behavior can be 
ruled out. On the other hand, in the absorbing phase per- 
sistence curves tend to a constant for large values of t, i.e. 
there is a non- vanishing asymptotic probability to persist 
indefinitely. This is a direct consequence of the fact that 
eventually the system reaches an absorbing configuration 
and the dynamics is arrested. Therefore, we conclude 
that unlike reversible systems where a power-law behav- 
ior holds for a wide range below the critical temperature, 
in systems with absorbing states the power-law behavior 
of the persistence is restricted to criticality. 

In what respect the global persistence, Figure 7 shows 
a double logarithmic plot of D g GB (t) versus t obtained 
at the second-order IPT for two different lattice sizes. As 
in the case of the local persistence, two time regimes ob- 
served: i) a short time regime (t < 300) where a power 
law behavior with a slope ~ 1.25±0.05, and ii) an asymp- 
totic regime (t > 300) with a power law behavior with 
exponent 9g GB ~ 2.5 ± 0.5. The onset of the long time 
regime depends on L as it is shows in figure 6 for L = 256 
and L = 512, however 9g GB appears to be size indepen- 
dent. As it was already discussed for the case of the local 
persistence, the onset of the crossover between the differ- 
ent time regimes can be linked with the results shown in 
figure 4. Consequently the asymptotic regime in the ZGB 
model starts when persistent sites are almost surrounded 
by inactive sites. 

Our result for the persistence exponents, namely 9i < 
9 g , is quite surprising at first sight, since from all pre- 
viously accumulated experience, mainly for equilibrium 
critical systems, it is known that 6i > 9 g , ||. It seems 
that irreversible critical systems may depart from this 
behavior. 

Assuming the dynamics of the order parameter to 
be Markovian, and using the known scaling relations, 
and values of the directed percolation universality class, 
one obtains that Eq. (1) can be written as [ p4|j33| 
9 g = 1 + d/(2z) w 1.565 (in d = 2). Our numerical 
value is far from this result, implying that the process 
is strongly non-Markovian. Oerding and van Wijland 
have shown p3], using field theoretical tools combined 
with the perturbative expansion method by Majumdar 
and Sire that up to first order in epsilon expansion, 
9 g = 2 + 0.059e + 9(e 2 ). This implies 9 > 2 (at least up to 
first order) as, in fact, is the case in our measurements. 

The physical interpretation for the fact that Qi < 9 g is 
qualitatively the following: at the critical point, in the 
activity density decay process it is more likely to fluctu- 
ate around the global average density (which corresponds 
to annihilating more (or less) activity than the average 
at that time), than invading absorbing regions in which 
persistent sites exist. This is very reasonable since in gen- 
eral, at the critical point and for large enough times, ac- 
tivity is restricted to some patches, which become smaller 
and smaller as time runs. It is therefore, more likely to 
have fluctuations in the density of these patches (con- 
trolled by the global persistence exponent) that to in- 



vade persistent sites (controlled by the local exponent) 
typically surrounded by absorbing regions at large times. 

The insert of figure 7 shows a log-log of Dg GB (t) . Here 
the events are short-lived making difficult the evalua- 
tion of the persistence exponent. However, considering 
larger error bars in this case, the exponent evaluated for 
Df- B (t) is consistent with the exponent previously cal- 
culated for Dg GB (t). So, within the discussed limitations 
for the evaluation of the exponents, we expect that both 
Dg GB (t) and Dg^ B {t) exhibit power-law behavior with 
the same exponent in agreement with the predictions of 
Oerding et al. ||. 

Since the ZGB model also exhibits a first-order IPT at 
Y2A — 0.52554, we have evaluated the persistence at the 
coexistence point. Figure 8 shows a log-log plot Df GB (t) 
versus t as obtained for lattices of size L = 512. Since 
the results are almost independent of L, finite size effects 
may be almost negligible. Figure 8 also shows that at 
coexistence the persistent sites are short-lived and conse- 
quently the persistence drops off abruptly. The data can 
be fitted by a power law with an exponent 9i ~ 3.0 ± 0.3. 
However, since the range of the fit is narrow and there 
seems to be a systematic curvature, we cannot disregard 
an exponential (or an stretched exponential) decay for 
longer times. Further clarification of this issue will have 
to wait for computationally expensive large-scale numer- 
ical simulations. It should be noticed that the operation 
of power-law (scale-invariance) in the dynamical critical 
properties of the first-order IPT of the ZGB model has 
recently been ruled out ||l[]. But given the lack of a 
general theoretical framework, we do not know whether 
or not algebraic decay of persistence functions should be 
expected at non-equilibrium first order transitions. 



V. CONCLUSIONS 

A numerical study of persistence in the ZGB dimer- 
monomer model for the catalyzed reaction A + — > 
AB2 is presented. It is found that for the second-order 
IPT the time dependence of both the local and the 
global persistence exhibit a crossover between a short- 
and a long-time regimes. A physical explanation for 
such a crossover has been provided. Persistence expo- 
nents are evaluated within the long-time regime giving 
9f GB ~ 1.50±0.01 for the local one, and 9 z g GB ~ 2.5±0.5 
for the global one. The former, which should correspond 
to the DP universality class in (2 + 1)— dimensions is very 
close to the value reported by Hinrichsen et al. HjJ for 
DP in (1 + 1) -dimensions, namely 9 BP ~ 1.50 ± 0.03. 
Simulations in the contact process (another model be- 
longing in the same universality class) confirm this in- 
distinguishability between the one and two-dimensional 
local persistence exponents with higher accuracy (see the 
striking parallelism of the curves in figure 5). This finding 
implies a theoretical and numerical puzzle: Is 9f p inde- 
pendent of the dimensionality for d < 2, or is it just that 
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the one and two-dimensional exponents incidentally take 
very similar numerical values? The likely possibility of 
them being equal is quite intriguing since the return prob- 
abilities in systems with absorbing states are expected to 
be rather different in different spatial dimensions. For 
the contact process, in dimensions larger than d — 2 we 
find, as expected, dimension dependent exponents. 

The fact that 8^ GB > 6f GB (in agreement with field 
theoretical predictions |24j] is in marked contrast to well 
established findings in the field of reversible transitions, 
where quite generically 9 g < B\. A physical explanation 
for such a discrepancy has been provided. 

We have also observed that power laws for persistence 
are obtained only at criticality, and the associated expo- 
nents do not depend upon the considered initial condi- 
tions, ad hoc with what is known for other well known 
reversible systems. 

Persistence probabilities have also been evaluated at 
the first-order IPT of the ZGB model. Due to the short- 
lived behavior of the persistent sites we obtained an expo- 
nent #f 3.0 ±0.3. However, further effort will be required 
in order to clarify if this apparent power-law behavior 
merely reflects a short-time regime followed by a cut-off 
or it is a real asymptotic effect. 

We expect our numerical findings to constitute a valu- 
able step for the development of a theoretical framework 
in the field of irreversible phase transitions, which is cer- 
tainly needed. 
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FIG. 1. Average coverage of B— species (8b), A— species, 
(6a) and the AB production rate (Rab), respectively, ob- 
tained under steady state operation, as a function of Ya 
for the ZGB model. Irreversible transitions of second- and 
first-order occurs at Yia ^ 0.387368 and Y2A — 0.52554, re- 
spectively 




FIG. 3. Log-log plot of Di(t) versus t, corresponding to 
the ZGB model at the second-order IPT and obtained for 
lattices of size L = 256. Results obtained starting simula- 
tions with different values of 9b (t = 0), as indicated in the 
figure, and averaging over 2.5 x 10 3 realizations. No significa- 
tive modification in the slope is observed upon changing the 
initial condition. 
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FIG. 2. Log-log plot of Di(t) versus t, corresponding to 
the ZGB model at the second-order IPT and obtained for lat- 
tices of size 0,L = 128 and A,L = 256 , respectively. The 
dashed (full) line shows the short- (long-)time regime and has 
slope 81 + I = 2.00 (8f GB + 1 = 2.5 ), respectively. Results are 
obtained starting simulations with 8s(t = 0) = 0.1 and aver- 
aging over 10 4 (2.5 x 10 3 ) realizations for L = 128 (L = 256), 
respectively. 
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FIG. 4. Linear-logarithmic plot of the average fraction of 
persistent sites with NN empty neighbors, for NN = 0, 1, 23, 
and 4 respectively, as a function of time t. Results are ob- 
tained using the same conditions as in figure 2. For times 
larger than 10 3 the fraction of empty sites neighboring per- 
sistent ones becomes very small; i.e. typically, activity takes 
place away from persistent sites, and local outbursts of activ- 
ity are required in order to reach persistent sites. 
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FIG. 5. Local persistence for one and two dimensional 
simulations of the contact process. In d = 1 (d = 2) the 
system size is L = 10 4 (256 * 256) respectively. Observe that 
both slopes are indistinguishable. The straight line has a slope 
6i = -1.5. 
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FIG. 6. Log-log plot of Di (t) versus t , corresponding to the 
ZGB model. Results obtained using lattices of size L = 256, 
with Quit = 0) = 0.10, and averaging over 2.5 x 10 3 real- 
izations. Different values of the parameter Ya were used, as 
indicated in the figure. The dashed (full) line shows the short- 
(long-)time regime and has slope 2.00 (0f GB + 1 = 2.5 ), and 
have been drawn to allow comparison with figure 2. See more 
details in the text. 
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FIG. 7. Log-log plot of Dg+(t) versus t corresponding to 
the ZGB model at the second-order IPT and obtained for lat- 
tices of side A,L = 256 and 0,L — 512, respectively. The 
dashed (full) line shows the short- (long-)time regime and 
has slope 1.25, (0g GB + 1 = 3.5), respectively. Results have 
been obtained starting the simulations with 6a {t = 0) = 0.10 
and averaging over 10 6 (2.5 x 10 s ) realizations for L = 256 
(L — 512), respectively. The inset shows a log-log plot of 
D g -(t) versus t, obtained for L — 512: the asymptotic slope 
coincides with that of the main plot. 




t 

FIG. 8. Log-log plot of Di(t) versus t obtained at the co- 
existence point (Y 2 a = 0.52554) for lattice sizes L — 256 and 
L = 512. The straight line has slope 9i + 1 = 4.00. Results 
have been obtained starting simulations with 9a (t = 0) =0.1 
and averaging over 10 4 (2.5 x 10 3 ) realizations for L = 256 
(L — 512), respectively. The upward bending of the last two 
points indicate that the value of the critical point we are con- 
sidering is slightly in the subcritical phase for the finite system 
sizes under analysis. 
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